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A new look at hadronic B decays 

Matthias Neubert ^ 

'^Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

We present a detailed study of non-leptonic two-body decays of B mesons based on a generalized factorization 
hypothesis. We discuss the structure of non-factorizable corrections and present arguments in favour of a simple 
phenomenological description of their effects. We discuss tests of the factorization hypothesis and show how 
it may be used to determine unknown decay constants. In particular, we obtain fn^ = (234 ± 25) MeV and 
fni = (271 ± 33) MeV. 



1. INTRODUCTION 

The weak decays of hadrons containing a heavy 
quark offer the most direct way to determine the 
weak mixing angles of the Cabibbo-Kobayashi- 
Maskawa matrix and to explore the physics of 
CP violation. However, an understanding of the 
connection between quark and hadron proper- 
ties is a necessary prerequisite for a quantitative 
theoretical description of these processes. The 
complexity of strong-interaction effects increases 
with the number of quarks appearing in the fi- 
nal state. Bound-state effects in leptonic decays 
can be lumped into a single parameter (the "de- 
cay constant" ) , while those in semileptonic decays 
are described by invariant form factors, depend- 
ing on the momentum transfer between the 
hadrons. Approximate symmetries of the strong 
interactions help us to constrain the properties of 
these form factors [Q. For non-leptonic decays, 
on the other hand, we are still lacking a com- 
prehensive understanding of strong-interaction ef- 
fects even in simple decay modes. The problem 
is exemplified in Fig. |l|, which shows multiple ex- 
changes of gluons between the quarks in the initial 
and final states. The intricate interplay between 
weak and strong forces has many surprising con- 
sequences, whose understanding is a challenge to 
theory. Examples are the A/ = ^ selection rule 
in K decays, the difference in the lifetimes of the 
charm mesons D"*" and and the difference of 
the lifetimes of the Af, and B particles. Although 
strong-interaction are more dramatic at low en- 
ergies, they are still hard to understand even in 
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Figure 1. Strong- interaction effects in a non- 
leptonic decay. 

B decays. 

At tree level in the Standard Model, non- 
leptonic weak decays are mediated by a single 
Ty-exchange diagram. When the external quarks 
have energies much below the electroweak scale, 
this process can be described by a local four- 
fermion interaction, which gets modified if hard 
giuon exchanges between the quarks are included. 
Their effects can be taken into account by using 
the renormalization group to evolve the effective 
interaction from the electroweak scale down to a 
scale fi of order mb, the mass of the decaying b 
quark. For the case of 6 — > cud transitions, e.g., 
the relevant part of the effective Hamiltonian is 

^ VcbV:,{c^{ii) {du)icb) + cifi) {cu)idb)} , (1) 

where (du) — ^7^(1 — 75)^ etc. are left-handed, 
colour-singlet quark currents. The Wilson co- 
efficients Ci(/i) are known to next-to-leading or- 
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der [^y^. At the scale /i — mt,, they have the 
values ci(mb) « 1.13 and C2(/i) « —0.29. These 
coefBcients take into account the short-distance 
corrections arising from the exchange of hard glu- 
ons. The effects of soft gluons (with virtualities 
below the scale /i) remain in the hadronic matrix 
elements of the local four-quark operators. A re- 
liable field-theoretic calculation of these matrix 
elements is the obstacle to a quantitative theory 
of hadronic weak decays. 

Using Fierz identities, the four-quark operators 
in the effective Hamiltonian may be rewritten in 
various forms. It is particularly convenient to re- 
arrange them in such a way that the flavour quan- 
tum numbers of one of the quark currents match 
those of one of the hadrons in the final state of 
the considered decay process. As an example, 
consider the decays B — > Dn. Omitting com- 
mon factors, the various decay amplitudes may 
be written as 

Ab«^d+.- = {c, + ^){D+7r-\{du){cb)\B°) 



Ar,o 



B<>^D"7r<> — C2 



+ C2 {D+TT-\^{dtaU){ctab)\B'>), 

^){D\\cu){db)\B') 



+ Ci {D\°\^{ctaU){dtab)\B'') , 
B-^DOtt- = Abo^d+tv- - V2Abo^£|0^o , (2) 



A 



where ta are the SU(3) colour matrices. The 
last relation follows from isospin symmetry of the 
strong interactions. The three classes of decays 
shown above are referred to as class- 1, class-2 and 
class-3, respectively [Q. 

2. FACTORIZATION HYPOTHESIS 

The above decay amplitudes contain the "fac- 
torizable contributions" 



{BD)-n: 



{Btx)D 



(3) 



The matrix elements in this equation are known 
in terms of the meson decay constants /^r and 
/£), and the transition form factors for the de- 
cays B ^ D and B ^ tt, respectively. Most 
of these quantities are accessible experimentally. 



Of course, the matrix elements appearing in 
also contain other, non-factorizable contri- 
butions. In general, we may define process- 
dependent hadronic parameters Si and es, which 
contain the non-factorizable corrections, in such 
a way that the decay amplitudes take the form 

AiB^'^D+n-) = aiJ^^BD)., 

A{B° ^D\") ^ a2T^B.)D, (4) 



with |-§ 
ci(m) 



ai 



0-2 



Nc 



(BD)tt 



■ C2(m) e 



C2(m) 



8 (P)^ 



Nc 



I 



_{Bt:)D 



(5) 



We stress that these expressions are exact; they 
are just a way to parametrize the relevant matrix 
elements of four-quark operators. The effective 
coefficients take into account all contributions 
to the matrix elements and are thus /z indepen- 
dent. The scale dependence of the Wilson coeffi- 
cients is exactly balanced by that of the hadronic 
parameters. 

Additional insight can be gained by combining 
these results with the 1/Nc expansion At the 
scale fi = 0{mb), the large- A^c counting rules of 
QCD imply ci = 1 + Oil/N^), ca = 0{1/Nc), 
= 0(1/7V2), and eg = 0(1/A^c). (For scales 
much lower than m^, the counting rules for the 
Wilson coefficients Ci (p) are spoiled by large loga- 
rithms.) Hence, we expect that <C 1, whereas 
contributions from eg can be larger. Using these 
results, we find from (0) 



ai - Climb) + Oil /N^), 
02 = C2(/i) + ciifi) + 4^''^^(/i) 
= C2(mh) + Cciimb) , 



0(1/^3) 
(6) 



where ( = 1/Nc + esirrib) is a process-dependent 
hadronic parameter of order 1 /Nc . It is impor- 
tant to stress that the naive choice oi — ci+C2/Nc 
and a2 = C2 + ci/Nc, which is often referred to as 
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"factorization hypothesis" , does not correspond 
to any consistent Umit of QCD; in particular, 
this is not a prediction of the l/N,. expansion. 
The more general expression for 02 given above 
was first introduced in Ref. jj]. Since the pa- 
rameter Eg is of order 1/Nc, the two contribu- 
tions to C are of the same magnitude, and hence 
C should be considered as an unknown dynami- 
cal parameter. As a general rule, we expect that 
non-factorizable corrections are small in class- 1 
transitions. In class-2 decays, on the other hand, 
the contribution proportional to eg is enhanced 
by the large value of the ratio ci/c2 — 0{Nc), 
and non-factorizable contributions can therefore 
be sizeable. 

The parameter eg obeys the renormalization- 
group equation 

f^^e^if^)^-^, (7) 

where the superscript P represents the depen- 
dence on the decay process. Let us assume that, 
for each process, there exists a "factorization 
scale" fif such that e^{fif ) = 0. (We will see 
later that this is indeed the case for two-body de- 
cays of B mesons.) Then (|^) implies that the 
phenomenological parameter ^ is given by 



Nc Stt fif ' 



(8) 



Based on the colour transparency argument of 
Bjorken we expect that fif scales with the en- 
ergy of the outgoing hadrons in a decay process. 
A fast-moving pair of quarks in a colour-singlet 
state acts as a colour dipole and decouples from 
soft gluons. Only hard gluons, with virtualities 
of order the energy of the outgoing particles, can 
rearrange the quarks and thus spoil factorization. 
On a qualitative level, the connection between the 
factorization scale and the energy release in the fi- 
nal state can be seen from Fig. H, where we show 
the ratio 02/01 as a function of as{fif). As we 
will see later, the value preferred by B ^ Dtt de- 
cays is positive and corresponds to a rather small 
coupling, indicating Hf — 0{Trn,) for these pro- 
cesses. On the other hand, D decays indicate a 
negative value of 02/01 , corresponding to a lower 
value of the factorization scale. This is in ac- 
cordance with the fact that in these processes the 
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Figure 2. The ratio O2/01 as a function of 
the running coupling constant evaluated at the 
factorization scale. The bands indicate the phe- 
nomenological values of 02/01 extracted from 
B Dtt and D — > Ktt decays. 

energy released to the final-state particles is much 
smaller. 

In two-body decays of B mesons, the energy 
release (per particle) in different processes differs 
by less than about 1 GeV. Combining this obser- 
vation with the empirical fact that in i? — > Dtt 
decays the scale fif is of order mi,, we find that 
the changes of the phenomenological parameter C 
in different decay channels are of order |^ 



AC 



few % . 



(9) 



Thus, on rather general grounds, we expect that 
in two-body B decays, to a good approximation, 
C is process independent. We shall refer to this 
assumption, together with the prediction for the 
parameters oi and 02 given in (^), as generalized 
factorization hypothesis. 

3. TESTS AND PREDICTIONS 

Adopting the theoretical framework described 
above, a prediction of hadronic decay amplitudes 
needs as input information meson decay constants 
and form factors. The decay constants of many 
light mesons are known with high accuracy from 
hadronic r decays, and from the electromagnetic 
decays of vector mesons [Q. The meson form 
factors of quark currents can, to some extend. 
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be extracted from seniileptonic decays. Whereas 
the various B D^*'' transition form factors are 
weh-known by combining data from semileptonic 
decays with heavy-quark symmetry relations |^], 
the available information about form factors de- 
scribing the decays of B mesons into light mesons 
is largely model dependent. As a consequence, 
for all class- 1 decays considered below the factor- 
ized decay amplitudes can be predicted without 
any model assumptions; however, the theoretical 
predictions for class-2 amplitudes involve larger 
theoretical uncertainties. In this work, we will 
adopt the NRSX model [0 to calculate the me- 
son form factors. Because of lack of information, 
we shall neglect final-state interactions between 
the produced hadrons in the final state. In ener- 
getic two-body decays of B mesons, these effects 
are expected to be small ||7|,pl|. 

Our predictions for the branching ratios of 
some of the dominant non-lcptonic two-body de- 
cays of B mesons are given in Table ||. The 
QCD coefficients ai and 02, as well as the (poorly 
known) decay constants of charm mesons, have 
been left as parameters in the theoretical expres- 
sions. For comparison, we show the world average 
experimental results for the branching ratios, as 
recently compiled in the review article in Ref . [l2| . 

3.1. Extractions of ai 

There are various ways in which to test the 
generalized factorization hypothesis. The most 
direct one relies on the close relationship between 
semileptonic and factorizcd hadronic decay am- 
plitudes. Comparing the non-leptonic decay rates 
with the corresponding differential semileptonic 
decay rates evaluated at the same value of pro- 
vides a direct test of the factorization hypothesis 
il. We have 



dT{BO^D(*Hi>)/dq^ 



(10) 



where h denotes a light meson {h = t: ov g), and 
fh its decay constant. To determine this ratio 
experimentally, one needs the values of the differ- 
ential semileptonic decay rate at various values of 
. They have been determined for B — > D*£D 



decays in Ref. fl^ , using a fit to experimental 
data. A comparison of the theoretical prediction 
with experimental data yields ai = l.llztO.lO for 
h = TT, and ai = 1.09 ± 0.13 for h — g. 

A more precise determination of the parame- 
ter fli is obtained by comparing the theoretical 
predictions for the decays B^ D^*'>~^h~ (with 
h = IT 01 g) directly with the data. For these pro- 
cesses, the transition form factors and decay con- 
stants are well known experimentally. Therefore, 
the theoretical uncertainties are minimal. From 
a fit to the data, we obtain ai = 1.08 ± 0.04. 

An interesting alternative extraction of ai 
is obtained from the class of decays S'^ 

, which involve the same 
transition form factors, but have a quite differ- 
ent kinematics since there is less energy release 
to the final-state particles. Unfortunately, the de- 
cay constants of Ds mesons arc not well known 
experimentally. From a fit to the data, we find 
ai = 1.10 ±0.07 ±0.17, where the last error takes 
into account the uncertainty in the decay con- 
stants. 

The different determinations of the parame- 
ter oi agree well with each other and yield a 
result that confirms the theoretical expectation 
that Oi ~ Ci(mb) « 1.13. Within the experi- 
mental errors, there is no evidence for a process 
dependence of ai. 

3.2. Extractions of 02 

A determination of the parameter 02, as well 
as of the relative sign between 02 and ai , is ob- 
tained by comparing the theoretical predictions 
for the decays B- D'-*'^°h- (with h = TT or 
g) with the data. Since the contributions to the 
decay amplitudes proportional to a2 involve the 
B ^ h form factors, they cannot be predicted 
in a model-independent way. Using the NRSX 
model for these form factors and assigning a con- 
servative error, we find 02/01 = 0.21±0.05±0.04, 
where the second error accounts for the model de- 
pendence. Combining this with the value for ai 
determined above gives 02 = 0.23 ± 0.05 ± 0.04. 
The fact that the ratio 02/01 is positive in B de- 
cays implies that class- 1 and class-2 decay ampli- 
tudes interfere constructively. This is in contrast 
with the situation encountered in charm decays, 
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Table 1 

Theoretical predictions for the branching ratios (in %) of non-leptonic B decays. In the third column, 
the factors containing not well known decay constants are suppressed. 



B° Modes 


NRSX Model 


ai = 1.08 
02 = 0.21 


Experiment 


Class- 1 


D+ir- 

D+Q- 

D*+-K- 

D*+Q- 

D+Dj 

D+D*- 

D*+D- 

D*+D*- 


0.257 a'i 
0.643 af 
0.247 al 
0.727 al 

0.879 al (/z5j240)2 
0.817 al (/z3./275)2 
0.597 al (/d]/240)2 
2.097 al (/d:/275)2 


0.30 
0.75 
0.29 
0.85 
1.03 
0.95 
0.70 
2.45 


0.31 ±0.04 ±0.02 
0.84 ±0.16 ±0.07 
0.28 ±0.04 ±0.01 
0.73 ±0.15 ±0.03 
0.74 ±0.22 ±0.18 
1.14 ±0.42 ±0.28 
0.94 ±0.24 ±0.23 
2.00 ± 0.54 ± 0.49 


Class-2 


K'>J/iP 

K*°J/iP 
K*°ip{2S) 


2.262 a'l 
1.051 al 
3.645 aj 
1.939 al 


0.10 
0.05 
0.16 
0.09 


0.075 ±0.021 

< 0.08 
0.153 ±0.028 
0.151 ±0.091 


B- Modes 


NRSX Model 


ai = 1.08 
a2 = 0.21 


Experiment 


Class- 1 


D°D*- 
D*°D- 
D*°D*- 


0.938 af ifDj240y 
0.873 0.1 (/z5./275)2 
0.639 al (/dJ240)2 
2.235 a?(/c./275)2 


1.09 
1.02 
0.75 
2.61 


1.36 ±0.28 ±0.33 
0.94 ±0.31 ±0.23 
1.18 ±0.36 ±0.29 
2.70 ±0.81 ±0.66 


Class-2 


K-J/ip 
K-^pi2S) 

K*-J/iP 
K*-ip{2S) 


2.411 al 
1.122 al 

3.886 al 
2.070 ai 


0.11 
0.05 
0.17 
0.09 


0.102 ± 0.014 
0.070 ± 0.024 
0.174 ±0.047 
< 0.30 


Class-3 


D°e- 

D*\- 


0.274 [ai + 1.127a2 (/d/200)]'' 
0.686 [ai + 0.587 aa (/d /200)]^ 
0.264 [ai + 1.361 az (/z3./230)]^ 
0.775 [al + 0.661 al (/d-/230)2 


0.48 
0.99 
0.49 
1.19 


0.50 ±0.05 ±0.02 
1.37 ±0.18 ±0.05 
0.52 ±0.08 ±0.02 
1.51 ±0.30 ±0.02 



+ 1.518 aiaa (/d-/230)] 
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where a similar analysis yields oi = 1.10 ± 0.05 
and a2 — — 0.49 ± 0.04 |10|, indicating a strong 
destructive interference. Since most D decays are 
(quasi) two-body transitions, this effect is respon- 
sible for the observed lifetime difference between 
D+ and D° mesons: t{D+) > t{D°). In B de- 
cays, on the other hand, the majority of tran- 
sitions proceeds into multi-body final states, and 
moreover there are many B~ decays (such involv- 
ing two charm quarks in the final state) where 
no interference can occur. The relevant scale 
for multi-body decay modes may be significantly 
lower than mt,, leading to destructive interference 
(see Fig. ||). Therefore, the observed constructive 
interference in the two-body modes is not in con- 
flict with the fact that t{B-) > t{B°). 

An alternative determination of the magnitude 
(but not the sign) of 02 can be obtained from the 
class of decays B K'^*'>tp^'\ which are charac- 
terized by a quite different decay kinematics. We 
find 02 = 0.21 ± 0.01 ± 0.04. A comparison of 
this result with the value of 02 determined above 
provides an interesting test of our theoretical pre- 
diction that even in decay modes with different 
energy release the process dependence of 02 is ex- 
pected to be mild. Within errors, there is indeed 
no evidence for any process dependence. Hence, 
the data fully support the generalized factoriza- 
tion hypothesis. From the result for 02 we may 
then extract the value of the phenomenological 
parameter ( or, equivalently, of the colour-octet 
matrix element eg. We obtain 

C = 0.45 ±0.05, £8(^6) = 0.12 ±0.05. (11) 

It would be most interesting to derive these re- 
sults from a rigorous, field-theoretical evaluation 
of the four-quark operator matrix elements. 

3.3. Determination of decay constants 

As an application, we shall employ the gener- 
alized factorization hypothesis to obtain rather 
precise values for the decay constants of the Ds 
and D* mesons. To this end, we derive from Ta- 
ble 1^ the theoretical predictions for the following 
ratios of decay rates: 



riB° 



D*+D: 



r(BO - 


D* + TT-) 


r(so - 


D+D*-) 


r(B" 


-^D+g-) 


r{B° - 


D*+D*-) 


r(BO - 


D*+Q-) 



0.72 



= 0.74 



1.68 



(12) 



r(B" 



D+D: 



T{B^ D+TT- 



1.01 



Theoretically, these predictions are rather clean 
for the following reasons: first, all decays involve 
class-1 transitions, so that deviations from fac- 
torization are probably very small; secondly, the 
parameter ai cancels in the ratios; thirdly, the 
two processes in each ratio have a similar kine- 
matics, so that the corresponding decay rates are 
sensitive to the same form factors, however eval- 
uated at different values; finally, we may hope 
that also some of the experimental systematic er- 
rors cancel in the ratios (however, we do not as- 
sume this in quoting errors below). Combining 
these predictions with the average experimental 
branching ratios ]T^ , 

B{B DD-) ^ (0.95 ± 0.24)% , 

B{B^DD*-) = (1.00 ±0.30)%, 

B{B^D*D-) ^ (1.03 ±0.27)%, 
B(B ^ D*D*-) = (2.26 ± 0.60)% , (13) 
we find the rather accurate values 
fjj^ = (234 ± 25) MeV , 

foi = (271 ± 33) MeV . (14) 

The result for is in excellent agreement with 
the value fo, = 241 ± 37 MeV extracted from 
the leptonic decay Dg fi^i' ffj]. The ratio 
fo'/fos — 1-16 ± 0.19, which cannot be deter- 
mined from leptonic decays, is in good agreement 
with theoretical expectations ||l5|,|l6|. Finally, we 
note that, assuming SU(3) breaking effects of or- 
der 10-20%, the established value of Jd^ implies 
Jd ~ 200 MeV, which is larger than most theo- 
retical predictions. 

4. SUMMARY 

Exclusive hadronic decays of B mesons are 
strongly influenced by the long-range QCD colour 
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forces. Theoretically, their description involves 
hadronic matrix elements of local four-quark op- 
erators, which are notoriously difficult to calcu- 
late. The factorization approximation is used to 
relate these matrix elements to products of cur- 
rent matrix elements. Conventionally, the factor- 
ized decay amplitudes depend on two phenomeno- 
logical parameters ai and a2, which are connected 
with the Wilson coefficients Ci (/i) appearing in the 
effective weak Hamiltonian. We have shown that 
this approach can be generalized in a natural way 
to include the dominant non-factorizable contri- 
butions to the decay amplitudes. (The situation 
is more complicated in B decays into two vector 
mesons, see Ref. |^].) In the generalized factoriza- 
tion scheme, the effective parameters ai and a2 
become process-dependent. However, using the 
large- A''c counting rules of QCD, we have argued 
that in energetic two-body decays of B mesons 
ai « ci(TOh) and a2 ~ C2(mb) -|- (ci{mi,), where 
( = 0{1/Nc) is a dynamical parameter. More- 
over, the colour transparency argument suggests 
that the process dependence of C is likely to be 
very mild, so that it can be taken to be a constant 
for a wide class of two-body decays. These the- 
oretical expectations are fully supported by the 
data. From a fit to the world average branch- 
ing ratios of two-body decay modes, we obtain 
ai « 1.08±0.04 and 02 ~ 0.21±0.05, correspond- 
ing to C ~ 0.45 ± 0.05. There is no evidence for a 
process dependence of these parameters; in par- 
ticular, the values obtained for 02 from the decays 
B and B- D(*)0/i-, where h = n 

or Q, are in good agreement with each other. 

We have discussed various tests of the general- 
ized factorization hypothesis by considering ratios 
of decay rates, and by comparing non-leptonic de- 
cay rates with semileptonic rates evaluated at the 
same value of q^. Within the present experimen- 
tal uncertainties, there are no indications for any 
deviations from the factorization scheme in which 
ai and 02 are treated as process-independent 
hadronic parameters. Accepting that this scheme 
provides a useful phenomenological concept, ex- 
clusive two-body decays of B mesons offer a 
unique opportunity to measure the decay con- 
stants of some light or charm mesons, such as 
the ai, Ds and D* 
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